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ABSTRACT 

About 30.5 days of nearly uninterrupted broadband photometry of the solar-type 
star k 1 Ceti, obtained with the MOST (Microvariability & Oscillations of STars) satel- 
lite, shows evidence for two large starspots with different rotation periods of 8.9 and 
approximately 9.3 days (Afi/fi ~ 4%). Ground based measurements in 2002 and 2003 
of Ca II H k, K emission reveal variations in chromospheric activity with a period of 
about 9.3 days. The data were obtained during the MOST commissioning phase. When 
the data are combined with historical observations, they indicate that the 9.3-day spot 
has been stable in its period for over 30 years. The photometry, with a sampling rate 
of approximately once per minute, was also used to search for acoustic (p-mode) oscil- 
lations in the star. We detect no clear evidence for p-modes in the k 1 Ceti photometry, 
with a noise level around 7-9 /umag at frequencies in the range 0.5 - 4 mHz (3-u 
detection limit of 21 - 27 //mag). There were no flares or planetary transits during the 
30.5 days of MOST monitoring with light amplitudes greater than 2 mmag (durations 
greater than 200 minutes) and 3 mmag (2 — 200 min durations). While this rules out 
any close-in planets of >0.5 Jupiter diameters with an orbital inclination close to 90°, 
the scatter in differential radial velocities permit a close giant planet in a more highly 
inclined orbit. 

Subject headings: stars: activity - stars: individual: Kappal Ceti - stars: late-type - 
stars: starspots - stars: rotation - stars: oscillations - stars: exoplanets 
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1. INTRODUCTION 

The production of solar flares, the evolution and migration of sunspots, and the very origin 
of the Sun's magnetic field are all believed to be associated with differential rotation beneath 
and at the solar surface. To test models of stellar dynamos, measurements of differential rotation 
in other stars are necessary for correlation with other parameters like magnetic variability and 
chromospheric activity. 

Surface differential rotation (SDR) can be easily observed in the Sun - with detailed obser- 
vations going back to Carrington (1860) - and subsurface rotation has been inferred from helio- 
seismic data of the Sun's 5-minute p-mode oscillations (e.g., Thompson et al. (1995)). However, it 
is a challenge to observe SDR directly at the surfaces of other stars. With the notable exception 
of EK Dra, where two spot modulations differing in period by 5.5% are simultaneously visible 
(Messina k, Guinan 2003), evidence for SDR in most solar- type stars with prominent starspots is 
seen only by comparing data over many years. Light variations ascribed to rotational modulation 
sometimes change in period from epoch to epoch, and the most natural explanation is that the 
dominant spot systems appear at different latitudes and move with different rotation periods (e.g., 
Hall & Busby (1990); Henry et al. (1995)). Doppler Imaging of some solar- type stars from epoch to 
epoch has also revealed evidence for SDR (see reviews by Collier Cameron (2002) and Strassmeier 
(2004)) and in at least one case, the rapidly rotating K dwarf AB Doradus, SDR was observed in 
the Doppler maps in only a few consecutive rotations (Collier Cameron et al. (2002)). So far, there 
have not yet been measurements of solar-type acoustic oscillations in any star other than the Sun 
with sufficient frequency resolution to explore its interior rotation. 

The MOST (Microvariability & Oscillations of STars) satellite (Walker et al. 2003a) was pointed 
at k 1 Ceti as part of the commissioning phase of the mission, to obtain useful science during 
engineering tests and debugging of the MOST systems. Our space-based photometry has sufficient 
time coverage and precision to reveal SDR which is easily evident even in only three equatorial 
rotations of the star. The data also set the first meaningful limits on the amplitudes of p-mode 
oscillations in this star. 

k 1 Ceti was chosen because of (1) its brightness and location in the sky at the time of these 
tests, (2) previous indications of spot migration and hyperactivity, and (3) its G5V spectral type, 
which made it a possible candidate for solar-type oscillations; we describe the star more fully in 
Section 2. The MOST observations are discussed in Section 3. The MOST-orbit binned data 
produce a light curve which reveals the flux rotational modulation of a young active Sun-like star 
in unprecedented completeness and precision (Section 4). We are able to relate the photometric 
modulations of spot visibility to periodic rotational variations of chromospheric activity seen in 
high-resolution spectroscopic observations of the Ca II K emission in 2002 and 2003 (Section 5), 



1 Based on data obtained with the MOST satellite: a Canadian Space Agency mission, jointly operated by Dynacon 
Inc., the University of Toronto Institute of Aerospace Studies and the University of British Columbia. 
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the latter obtained just before the MOST observations; the same data provided a new accurate 
determination of V sini for the star. The full temporal resolution MOST photometry is used to set 
meaningful limits on the oscillation amplitudes in this star in Section 6. 

2. THE STAR k 1 Ceti 

k 1 Cet (HD 20630, HIP 15457, HR 996; V = 4.83, B — V = 0.68) is a nearby (9.16 ± 0.06 pc, 
ESA (1997)) G5V dwarf. Its variability with a period of 9.09 days was detected by the Hipparcos 
mission (ESA 1997). Since then, several studies were aimed at reconciling apparent changes in the 
period, changes which can be explained by different latitudes of spot formation in different years. 

As a MK spectral classification standard (Gray et al. 2001), k 1 Ceti is one of the most fre- 
quently observed stars. Although sometimes considered a "Very Strong Lined" late-type dwarf, its 
metallicity may be only slightly higher than solar, [Fe/H] = +0.05 ± 0.05 (for full references, see 
Heiter k, Luck (2003)). As far as it has been possible to establish, it is a single star and does not 
possess any large planets (Halbwachs et al. 2003). Its radial velocity of +18.9 km s _1 combined 
with the Hipparcos proper motions leads to a rather moderate spatial velocity relative to the Sun, 
suggesting Young Disk population membership (the two available estimates disagree slightly: Gai- 
dos et al. (2000), U = -12.7, V = +7.1 and W = +2.6; Montes et al. (2001) U = -22.4, V = -4.3 
and W = -5.3). 

Giidel et al. (1997) estimated an age of 750 Myr from the relatively rapid rotation of 9.2 days 
seen in the spot modulation and suggested that the star is a likely member of the Hyades moving 
group. However, Montes et al. (2001) considered seven moving groups of nearby young stars, but 
were unable to associate k 1 Ceti with any of them. 

The young age of the star was the reason for Dorren & Guinan (1994) to include k 1 Ceti in 
the The Sun In Time project which attempts to characterize young solar-type stars in terms of 
temporal changes taking place particularly at the epochs before terrestrial life formation (Guinan 
et al. 2003). In the group of 6 such stars, with the youngest being 130 Myr old, k 1 Ceti is one of 
the most advanced ones with an estimated age in the Dorren & Guinan (1994) study of about 650 
Myr. The difference in the age of 100 Myr versus 750 Myr of Giidel et al. (1997) can be traced to 
the current uncertainty of estimates based on the rotation rate. (Note: The 8.9 day rotation period 
found in this paper may be taken as an indication of an even younger age). 

As is observed for stars with activity somewhat moderated with age, the star shows an ac- 
tivity cycle. Baliunas et al. (1995) monitored the narrow-band Ca II H & K chromospheric fluxes 
photoelectrically from 1967 to 1991 expressing them in terms of the Shk index. They found a 
rotational period of 9.4 ± 0.1 days (Baliunas et al. 1983) with a chromospheric activity cycle of 5.6 
years (Baliunas et al. 1995); the quality of the latter determination was described as "fair" and 
longer term trends were noted. Using broad band photometry, Messina & Guinan (2002) observed 
a photometric activity cycle of of 5.9 ± 0.2 years between 1990 and 1999, with a rotational period 
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of 9.214 days. 

On 24 January 1986, Robinson & Bopp (1988) caught the signature of a massive flare in 
He I D 3 which Schaefer et al. (2000) subsequently estimated had an energy ~ 10 35 ergs. Schaefer 
et al. (2000) speculated that such massive flares on a single solar-type star might be triggered by 
the magnetic field of a close giant planet such as those in 51 Peg- type systems. Neither Walker et 
al. (1995) nor Cumming et al. (1999) detected any long-term periodic variability over many years 
in their precise radial velocity (PRV) measurements, so there exist no dynamical perturbations 
consistent with a close-in giant planet (unless the orbit is highly inclined). On the other hand, 
Walker et al. (1995) did find (their Figure 2) a rapid RV change of +80 m s _1 in 1988.5 which 
was accompanied by an equally rapid increase in chromospheric activity in the Ca II A8662 line. 
While the RV change could be modeled by the close approach of a giant planet in a highly elliptical 
orbit, the tight correlation of the RV with chromospheric activity (both positive and negative RV 
excursions) pointed to a change intrinsic to the star. In Section 5.4 we point out that, within the 
scatter of existing RV measurements, k 1 Ceti could harbor in a close, high inclination orbit. 



3. THE MOST PHOTOMETRY 

3.1. Photometry of k 1 Ceti 

The photometric observations discussed in this paper were obtained by the MOST (Microvari- 
ability & Oscillations of STars) microsatellite (Walker et al. 2003a) between November 5 and De- 
cember 5, 2003, during the commissioning phases of the mission (constraints and limitations of these 
phases are described in the next Section 3.2). MOST is Canada's first orbiting space telescope, 
launched in June 2003 to study (1) acoustic oscillations in solar-type and magnetic peculiar stars, 
(2) rapid variability in hot massive stars; and (3) reflected light from close-in giant exoplanets. The 
MOST instrument is an optical telescope (aperture 15 cm) feeding a CCD photometer through a 
single broadband optical filter (350 - 700 nm). Thanks to its polar Sun-synchronous orbit (altitude 
820 km), MOST can monitor stars in a zodiacal band of sky about 54° wide for up to 8 weeks 
without interruption. Starlight from Primary Science Targets, such as k 1 Ceti, i.e. brighter than 
V ~ 6, is projected onto the CCD as a fixed image of the telescope pupil covering about 1500 
pixels for high photometric stability and insensitivity to detector flatfield and radiation effects on 
individual pixels. As has already been demonstrated for Procyon (Matthews et al. 2004), MOST 
is capable of reaching in that mode a photometric precision of about ±1 part per million (ppm) to 
search for oscillations in very bright stars. 

Despite limitations during commissioning (see the next section), the MOST data of k 1 Ceti 
achieved a duty cycle of about 96% over a time span of 30.5 days (with only three gaps of a few 
hours each), with 40-second exposures obtained approximately once per minute. At the time these 
data were obtained, this was the longest uninterrupted observational coverage of any star other 
than the Sun. (MOST has since exceeded this duty cycle and time coverage on other targets, in- 
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cluding Procyon: Matthews et al. (2004).) MOST was designed to be highly stable photometrically 
over periods of hours and self-calibrating (see the discussion in Walker et al. (2003a)), consequently 
the measurements are non-differential. For some Primary Science Targets, fainter Secondary Sci- 
ence Targets in the adjacent direct-focus field may allow us to search for systematic drifts in the 
photometry, but no such stars were observed with a sufficient S/N for k 1 Ceti. 

The MOST photometry is performed in a single spectral band, with a filter about twice as 
wide as the Johnson V filter and slightly redward of it (Walker et al. 2003a) , without reference to 
comparison stars or flux standards. Thus, there is no direct connection to any photometric system 
and hence, no absolute calibration of the MOST photometry. Because of this, the original data 
(expressed in Analog-to-Digital intensity units or ADU) were normalized to unity at the maximum 
observed flux of the star, occurring at JD 2,452,953.0. All variation amplitudes in this paper are 
also expressed in these relative units. 

To investigate rotational modulation in k 1 Ceti, the 40-second exposures were binned according 
to the 101.4-min orbital period of the MOST satellite, increasing the signal-to-noise per data point 
and minimizing the influence of orbital variations in the photometric background due to stray 
Earthlight. The electronic-version Table 1 lists the MOST orbit-averaged photometry, and those 
points are plotted in the upper panel of Figure 1. 

The formal mean standard error of a single data point is 0.77 mmag (770 parts per million or 
ppm), which is about twice the error calculated from the detected photon shot-noise. The simple 
average of 101 individual 40 sec integrations would give an error of ~100 ppm, which is not far from 
what is observed as the actual scatter of about 200 ppm - which includes the intrinsic variability of 
the star. Some of the photometric noise is certainly intrinsic to the star (e.g., granulation variations, 
given the star's activity) and some is undoubtedly due to a non-Poisson instrumental noise resulting 
from the significant pointing scatter and the removal of the (noisy) stray light background signal 
(see the next Section). 

The sky background removed from the data shown in Figure 1 varied from about 21 to 25% 
of the maximum flux measured from the star (on JD 2,452,953); the background variations are 
plotted in Figure 2. The background light is due primarily to light entering the instrument from the 
illuminated limb of the Earth which is exposed to the face of the satellite with the telescope entrance 
aperture. During commissioning, the MOST team was still investigating the background and 
exploring ways to mitigate it for future targets (e.g., choosing the optimum "roll" orientation for the 
satellite). To correct for the sky background, which is modulated with the satellite orbital period, 
but varies slightly from orbit to orbit, a running averaged background phased with the orbital 
period was subtracted from the data. During the k 1 Ceti observing run, the Moon approached 
within about 14° of the MOST target field, but we observe no correlation of the sky background 
with the angular separation of the Moon and the target, or with the lunar phase (see the lower 
panels of Figure 2). The main results of this paper are not particularly sensitive to any residual 
background variations which may be present. Although the sky background variability is probably 
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the fundamental limit to the long-term precision of the K Ceti photometry, we note that the 8.9- 
day period light photometric cycle was highly stable in the final light curve (Section 4.2), indicating 
that the background was correctly accounted for in the MOST data over the duration of the k 1 Ceti 
observing run. 

3.2. Nature of the Commissioning Science observations 

A brief description of the "commissioning" observations is given below as a record and expla- 
nation of certain limitations of these observations compared to the normal scientific operations of 
MOST which were initiated in January 2004, not long after the k 1 Ceti data were obtained. 

The precise photometry that MOST was designed to obtain is based on integrating the signal 
in a large image of the telescope pupil (covering about 1500 pixels of the CCD) projected by a Fabry 
microlens and field stop enclosing the incoming stellar beam. In normal operations, a resolved sub- 
raster of the CCD with that pupil image is downloaded for every science exposure, along with Fabry 
images of the adjacent sky backgrounds, for processing on the ground. This is known as Science 
Data Stream 2 or SDS2. As a data backup, we also perform the pupil image signal integrations 
on board the satellite and store that information there for several days in case of extended loss of 
contact with our ground station network. These data are part of Science Data Stream 1 (SDS1). 

During commissioning, there was an unrecognized problem with a small fraction of the onboard 
memory which led to computer crashes occurring roughly every 1 — 2 days if both SDS2 and 
engineering telemetry were being sent to Earth. This resulted in a fraction of about 21% data 
sent as the fully resolved SDS2 files; the rest of our k 1 Ceti photometry is based on SDS1 data. 
The lack of spatially resolved Fabry image data restricts our options of sky background and stray 
light removal. However, our later experience (see Matthews et al. (2004), for the example of 
Procyon) shows that the quality of SDS1 photometry is comparable to SDS2, although the latter 
is preferred. Even with our routine of downloading only SDS1 data, the onboard memory problem 
led to occasional computer crashes, which interrupted the data flow for one or more satellite orbits 
at a time, and accounts for most of the gaps in the light curve of Figure 1. The small percentage 
of corrupt memory locations has since been identified, and we avoid storing data there, so this 
problem no longer affects the MOST mission. 

At the commissioning phase of the mission, we had not yet optimized the spacecraft pointing. 
The mean pointing errors were 4.5 arcsec rms in x and 6.2 arcsec in y, so this image wander across 
the Fabry microlens (with a field stop radius of 30 arcsec) led to larger photometric errors than 
were later obtained during normal scientific operations. Also, this limited our ability to center 
the target within the field stop. Jitter of the star within the field stop increases noise. These are 
no longer serious issues for MOST. Upgrades to the Attitude Control System software have led 
to steady improvements in pointing. Observations of Procyon in January - February 2004 were 
made with pointing errors of 1.3 and 3.1 arcsec rms in x and y respectively; the current pointing 
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performance has been improved to about 1.0 arcsec rms in both axes. 

4. SPOT VARIABILITY AND EVIDENCE FOR SURFACE DIFFERENTIAL 

ROTATION 

4.1. Photometry of k 1 Ceti before MOST 

The current explanation of the variability observed in k 1 Ceti is that stellar spots cause light 
changes at the level of 0.02 to 0.04 magnitude. The changes have been observed by several ground- 
based observers to have characteristic time scales of about 9 days. The most extensive studies were 
done by Gaidos et al. (2000) and Messina & Guinan (2003). From 20 years of literature data, Gaidos 
et al. (2000) found several periods ranging between 9.14 and 9.46 days, each with uncertainty at the 
level of 0.03 - 0.05 days. Messina & Guinan (2003) conducted a careful study over 6 seasons and 
determined individual periods ranging between 9.045 and 9.406 days, with typical uncertainties of 
0.02 - 0.05 days. The simplest explanation for such period changes is that spots are forming at, or 
migrating to different stellar latitudes, rotating with different angular velocities. 

By relating the spot modulation periods to the photometric cycle, Messina & Guinan (2003) 
suspected that k 1 Ceti shows an "anti-solar" pattern in which the period of the rotation modulation 
tends to increase steadily during the activity cycle. 

4.2. Interpretation of the MOST photometry 

The photometric data discussed here consist of the orbit-averaged observations (binned by 
intervals of 101.4 min) corrected for sky background, as presented in Table 1 and the upper panel 
of Figure 1. A double- wave variation is clearly evident, consisting of deep minima with a depth of 
about 0.04 of the maximum light, interspersed with moderately shallow depressions with amplitudes 
of about 0.005 — 0.01. The shallow minima progressively move relative to the deeper ones as if two 
close periodicities were involved. 

The photometry only covers about three cycles of the dominant variation, so Fourier analysis 
is relatively ineffective at accurately identifying the periodic content of the variations. We formally 
recover two significant peaks at frequencies of 0.1146 and 0.2170 day -1 , corresponding to periods 
of 8.72 and 9.22 days (actually 4.61 = 1/2 x 9.22). Because of the short duration of the observing 
run, these periods cannot be established to better than about ±0.1 day (Figure 3). Instead of using 
the formal values of the derived periods, we use two periods, 8.9 and 9.3 days, which are derived 
by a simple folding of the data in a sequential decomposition of the light curve, starting from the 
largest and most regular variation: 

1. The most obvious is a periodic pattern reminiscent of an eclipsing light curve or a spot curve, 
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the latter due to a spot or a group of spots with probably well denned edges covering a 
moderately large fraction of the stellar surface. This pattern has an amplitude reaching 0.04 
of the maximum observed flux and a period of 8.9 ±0.1 days. An average of three cycles gives 
a relatively stable light curve, as shown in the phase diagram of Figure 4. 

2. When the 8.9 day pattern is removed, another periodic variation becomes visible. This is 
shown in the middle panel of Figure 1. This pattern is less regular than the 8.9-day variation, 
with the minima becoming progressively deeper from about 0.005 to 0.01 over the three 
complete cycles. (There is evidence that the next cycle, which is incomplete, is shallower 
than the preceding one.). The three times of minima give an approximate period of 9.3 days, 
where the uncertainty is harder to estimate due to the variable shapes of the dips. 

3. The middle panel of Figure 1 suggests that either the 9.3 day pattern evolved over time or 
there existed residual variability from some other source, especially at the beginning of the 
MOST run. The variability was most probably intrinsic to the star, judging by the progression 
of the depths of the 9.3 day pattern, and the excellent agreement between the phased minima 
of the 8.9 day pattern. We observe no correlation in the sky background variations, shown 
in Figure 2, with the deviations from regular 9.3 day periodic pattern, shown in the middle 
panel of Figure 1; this may be taken as indication that the sky background variations were 
not the cause of these irregularities. 

4.3. Spot modeling 

The observed periodic variations can be explained by spots formed at different stellar latitudes 
and carried by the rotation of the star with different angular velocities. With single-band photom- 
etry, little can be said about the type of spots, if they are hotter or cooler than the surrounding 
photosphere. However, two of the authors did develop two simple, entirely independent geometric 
spot models to try reproduce the observed light curve. The models differed in approach in that 
the first used the minimum number of parameters and ended up describing only the larger spot 
which apparently did not evolve during the MOST observations; only a very simplified description 
could be derived for the smaller spot. The second model utilized a full set of 12 parameters re- 
quired to describe two spots and the geometry of their visibility on the star. The second model 
was under-constrained, but its basic results fully confirmed the first model in its description of the 
larger spot. 

In the simplest terms: The larger of the two spots, with a rotation period of 8.9 days and a 
relative flux amplitude of 0.039, if entirely black and circular, would have a radius of 19% of the 
stellar radius or an angular radius from the star center of 11°. The mean 8.9-day light curve is 
shown in Figure 4. Since the light variations take place over slightly more than half of the rotation 
and the "flat top" is shorter than half of the rotation, the spot appears to be situated above the 
equator. The smaller spot, causing the 9.3-day (or slightly longer) periodicity, also if black and 
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circular, must be about 12% of the stellar radius (angular radius 7°). 

Any model must face the limitation that the parameter without any constraint is the unper- 
turbed level of the light curve, Iq; we arbitrarily assumed it to be at the highest observed level (see 
Figure 1), but some small spots may have been present on the surface at all times. While the spot 
longitude (A) can be removed from the parameter list by an appropriate adjustment of the phases, 
even the simplest model involves several adjustable parameters: the rotation axis inclination (i), 
the limb darkening coefficient (u), the latitude (/?), the size or the angular radius as seen from 
the center of the star (r), as well as the darkness (i) of the spot. The darkness (0 < t < 1) of 
the spot cannot be estimated from the single-color data, so that the minimum size of the spot 
can be estimated assuming a black spot, t = 0. Even such a highly simplified model, with several 
parameters fixed contains at least three fully adjustable geometrical parameters: i, (3, r. The best 
model of the larger spot fits the primary variation to better than 0.001 (Figure 4). The relevant 
parameters for the star and the larger spot are (the symbol = means a fixed parameter): 

P = 8.9 ± 0.1 day, i = 70° ± 4°, (3 = +40° ± 7°, r = 11.1° ± 0.6°, A = 0°, t = 0, u = 0.8 

The 8.9-day spot is relatively large, and may be larger if not entirely black, so that the assumption 
of it being circular - especially in view of the clear evidence of the shearing differential rotation 
- is highly disputable. Therefore, this solution should be taken as indicative rather than of much 
physical significance. 

The second model assumed two independent spots with different rotation periods. It was 
specified by the rotational inclination i of the star and the limb darkening parameter u, with all 
parameters describing the two spots: the radii and flux contrasts, and the latitude, initial longitude, 
and the rotation period for each spot. Such a model required extensive tests of different combina- 
tions of parameters. While the results are inconclusive in details because of the large number (12) 
of weakly constrained parameters, they are fully consistent with the minimum-parameter model 
described above in the description of the first, 8.9-day spot. The second spot in the extended model 
appears to have a longer period, perhaps as long as 9.7 days, although uncertainty of the period is 
large and results from the obvious evolution of the second spot over time and a complex coupling 
with values of the remaining parameters. 

The models have been unable to define accurate latitudes of both spots, although the second, 
more complex model suggested a higher latitude for the smaller, 9.3 day spot, a situation which 
would correspond to the solar surface differential rotation pattern. However, there is no question 
from both independent analyses that the two spots have unambiguously different periods, support- 
ing the interpretation of SDR with AQ/£l ~ 4 %. This is a clear evidence of surface differential 
rotation on k 1 Ceti, and only the second after EK Dra (Messina & Guinan 2003) where two spots 
were also simultaneously visible. All other inferences were based on apparent period changes (Hall 
& Busby 1990; Henry et al. 1995; Gaidos et al. 2000; Messina & Guinan 2003), interpretation of 
subtle broadening effects in spectral lines, so far detectable only in the most rapidly rotating stars 
(Donati & Collier Cameron 1997; Donati et al. 1999) or Doppler Imaging from epoch to epoch (Col- 
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Her Cameron 2002); for general reviews of these techniques and their limitations, see Rice (2002), 
Gizon & Solanki (2004) and Strassmeier (2004). 

4.4. Searching for flare activity and transits 

The MOST orbit-averaged data were also searched for short-lived variations which might be 
associated with transits of a close-in giant planet or with white-light flares. The latter seemed more 
likely given the history of activity of k 1 Ceti. No candidates for transits or flares were found in the 
data. We can set useful limits on such events over the month that we observed the star. There were 
no flare events or transits with peak-to- valley variations > 2 mmag with durations longer than 200 
minutes and none with variations > 3 mmag with durations between 2 and 200 min. 

For planetary transits, this limits the size of any possible transiting bodies with orbital periods 
of about 30 days or less to approximately 0.045 — 0.055 of the stellar radius, about half the size 
of Jupiter. However, our spot modeling and measurements of V sin i indicate that the rotational 
inclination of the star is not close to 90° so it is not expected that the orbital plane of any planets 
in the k 1 Ceti system would be coincident with the line-of-sight. 

5. GROUND-BASED SPECTROSCOPY 

5.1. The CFHT spectroscopic observations 

Several of us (D.B., S.R., E.S., G.W.) monitored the Ca II H & K reversals of a number 
of 51 Peg systems (Shkolnik et al. 2003) as well as several active, single stars with the Gecko 
spectrograph on the Canada France Hawaii 3.6-m Telescope (CFHT). The intention was to detect 
chromospheric activity synchronized with the planetary orbits. As a by-product, we determined 
radial velocities (RVs) for these stars with orv <20 m s" 1 (Walker et al. 2003b). Here we present 
the various results for k 1 Ceti. 

The dates of the observations are listed in Table 2. They were acquired on four observing runs 
in 2002 and 2003 with the Gecko coude spectrograph, fiber-fed from the Cassegrain focus of the 
CFHT. Single spectra were recorded at a resolution of 110,000 spanning 60 A and centered at 3947 
A (see Figure 5). The signal-to-noise was ~ 500 per pixel in the continuum and 150 in the H & K 
cores. Full details of the data reduction are given in Shkolnik et al. (2003). 

5.2. Determination of Vsini 

k 1 Ceti shows moderate broadening of lines due to rotation at the rate of about 9 days. The 
broadening is detectable in our high-resolution spectroscopy. Combined with the rotation period 
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and the value of the radius, V s'mi can be used to set a constraint on the rotational axis inclination 
i which is one of free parameters in the spot models. The previous determinations of V sini were 
3.9 km s~ l (Fekel 1997), 4.3 km s" 1 (Benz & Mayor 1984), 5.6 km s" 1 (Soderblom et al. 1989) and 
3.8 km s" 1 (Messina et al. 2001). 

We determined the value of V sin i by analyzing the broadening of the photospheric lines in 
the region of Ca II 3933/3968 A in our CFHT spectra. We used the broadening function formalism 
(Rucinski 1999, 2002) of the deconvolution of photospheric spectra of rotating stars with spectra 
of sharp-line, slowly-rotating stars. The approach leads to removal of all common agents of line 
broadening, not only of the instrumental component, but also of the tangential turbulence and of 
all other components which are the same for the program and the template stars. For k 1 Ceti, we 
used r Ceti which is a slightly cooler dwarf, G8V, and rotates very slowly, Tasini = 0.6 km s _1 
(Fekel 1997). 

For determination of the broadening function, the spectra must be rectified and normalized. 
As shown in Figure 5, the rich photospheric spectrum in both re 1 Ceti and r Ceti is superimposed 
on strongly variable far wings of the Ca II absorption. We removed sections of the spectrum where 
the wings fall below 1/3 of the maximum flux and applied a simple, upper-envelope rectification 
between high points of the quasi-continuum. Figure 6 shows the broadening function for one of 
the CFHT spectra, obtained on August 26, 2003 (see Table 2). Instead of least-squares fits of the 
limb-darkened rotational profile (Gray 1992) to the broadening functions (with the limb darkening 
coefficient assumed u = 1.0, appropriate for the UV in late-type stars), we matched the half-widths 
at half maximum; such V sin i determinations are less sensitive to the residual uncertainties at the 
base of the broadening function. Our average Vsini = 4.64 ± 0.11 km s _1 is compatible with 
previous determinations. We note that Gaidos et al. (2000), using an almost identical value of 
Vsini = 4.5 km s _1 , interpreted it as implying 40° < i < 90°. The main spot in both model fits is 
located at moderate stellar latitude above the equator. If the period of 8.9 ±0.1 days were ascribed 
to the stellar equator, then the best fitting inclination of i = 60° ± 5° and the measured Vsini 
combine to give a stellar radius of R ~ O.95±O.lOi? , which is certainly reasonable for a G5V star. 



5.3. The K-line residuals 

The 7 A region centered on the K-line which was used in the analysis is shown for a single 
spectrum in Figure 7 where spectral intensity is normalized to unity at the wavelength limits. In 
Figure 8 the residuals from the average normalized spectra are shown for each night (where pairs of 
spectra on the same night have been averaged). Activity is confined entirely to the reversals, within 
about 0.94 A around the center of the emission core or within ±35.8 km s _1 . If this broadening 
is caused by the velocity field, it is some 7.7 times larger than Vsini. To visualize the changes 
in the reversals, the integrated flux residuals shown in Figure 9 are expressed in A relative to the 
restricted continuum height shown in Figure 7. This height corresponds to about 1/3 of the local 
photospheric continuum around the Ca II lines so that the integrated residuals are proportional, 
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but not identical, to the emission equivalent width. They are listed in column 4 of Table 2. 

The best fitting periodic variation for the K residuals in Table 2 over the four observing runs in 
2002 and 2003 has a period of 9.332 ±0.035 days. The residuals are plotted as a function of phase in 
Figure 9 with the arbitrary initial epoch of the minimum residual flux at To{JDh e i) = 2, 452, 485.117 
(the minimum flux epoch can be estimated to ±0.1 days). The data from 2002 are shown as open 
symbols and those from September 2003 are solid. Although the data from 2002 cover only 3 
rotations, with the data from 2003 the span in the rotations is extended to 44. 

Since the period is largely defined by the 2002 observations, Figure 9 implies that the rotation 
associated with the enhanced Ca II activity has remained quite stable. The 9.332 day period is 
entirely compatible with that found by Baliunas et al. (1983) (9.4 ± 0.1 day) suggesting that the 
enhanced activity has in fact been stable at the same latitude for at least 36 years (over 1400 
rotations). This period is also compatible with the 9.3 day period of the less prominent spot that 
we see in the MOST observations (see Section 4.2). Given this apparent stability, it is reasonable 
to extrapolate variations in the K-line residuals to the epoch of the MOST photometry which 
corresponds to only 7 additional rotations. The bottom panel of Figure 1 shows the K-line flux 
changes which can be related to the photometric variations due to the 9.3 day spot seen in the 
MOST photometry. The maxima of the two curves coincide. The interpretation of this coincidence 
is not obvious as the enhanced Ca II could be expected to accompany the dark spot, unless what 
we see as the 9.3 day photometric signal is actually a bright spot and the minima are due to 
disappearance of the bright plage region. 



5.4. The Radial Velocities 

Radial velocities were estimated from m 20 A of the photospheric spectrum (3942 — 3963 A, 
see Figure 5) with the "fxcor" routine in IRAF 5 , with Th / Ar comparison lines providing dispersion 
correction. More detail is given in Walker et al. (2003b). The first spectrum of the series in 2002 
acted as the template from which differential radial velocities (ARV) were measured for all the 
spectra. The ARV values are listed in Table 2. They are plotted against the 9.332 day rotational 
phase determined from the K-line residuals in Figure 9 where the open symbols correspond to data 
from 2002 and the solid points to September 2003. There is a remarkable consistency between the 
ARVs from both years suggesting that the underlying ARV technique is highly stable even without 
an imposed fiducial (such as iodine vapor). 

From the 2002 data we derive <jrv =21.8 m s _1 and 23.6 m s _1 when combined with 2003. 
These values are very similar to (Try =24.4 m s _1 , found over 11 years by Cumming et al. (1999). 



5 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association 
of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science 
Foundation. 
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For comparison, a planet which induces a reflex radial velocity variation < 50 m s -1 and is tidally 
synchronized to the star, would have Msini ~ 0.7AM j and a = 0.084 AU. Consequently, k 1 Ceti 
could still harbor a giant planet in a tidally locked orbit which we would not have detected. 

In 2003 the ARVs appear significantly different between the two nights, something which seems 
to be reinforced by the consistency within the pairs of ARVs. While a planetary perturbation cannot 
be ruled out by the 2002 data and other Precision Radial Velocity studies, the difference in 2003 
might be associated with the velocity field of the star itself. The increase of velocity with increasing 
K-line strength is consistent with the extreme event in 1988.5 seen by Walker et al. (1995). However, 
the data are insufficient to discuss the point further here. 

6. ABSENCE OF P-MODE OSCILLATIONS 

In Figure 10, we present the Fourier amplitude spectrum of the complete set of unbinned pho- 
tometry (at full temporal resolution of about one minute), in which the lower frequency modulations 
have been filtered from the data. Also, the known satellite orbital frequency and its harmonics have 
been removed from the data, to avoid any confusion with residual stray light variations. Although 
the spectrum is plotted only to a frequency of 4 mHz (period = 4.17 min), the noise is level out 
to the Nyquist frequency of about 8 mHz, and there are no significant periodic signals detected in 
that range. The white dashed line represents the theoretical Poisson noise for the data. 

Although there are a few peaks visible at frequencies below 1 mHz, there is no evidence for the 
pattern of nearly equal spacing associated with p-mode oscillations of low degree and high overtone 
(as are seen in the Sun). There are several isolated peaks in the range 0.1 - 1.0 mHz (periods of 
2.8 hr to 17 min) which may be real but are too few and widely separated to attempt a meaningful 
fit to a model. 

At a detection level of 3a, we exclude coherent oscillations larger than about 27 ppm around 
1 mHz, 26 ppm at 2 mHz, and 24 ppm at 3 mHz. (The last is in the frequency range of the solar 
oscillations.) Thus, we would not expect to detect p-mode oscillations in k 1 Ceti if they were at 
the level observed in photometry of integrated light of the Sun; i.e., about 4-10 ppm (Frohlich et 
al. 1997). 

However, knowledge of the nature of p-mode oscillations in other stars is very limited, and 
the first space-based photometry of Procyon Matthews et al. (2004) obtained by MOST yielded an 
unexpected null result despite a detection limit at least twice as good as in these commissioning 
data of re 1 Ceti. We are only beginning to explore the characteristics of p-mode oscillations in 
luminosity for stars beyond the Sun. The moderately young age of k 1 Ceti, its relatively rapid 
rotation compared to the Sun, and its later-than-solar spectral type may be contributing factors to 
the absence of large-amplitude (i.e.,> 20 ppm) p- modes in this star. 
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7. CONCLUSIONS 

We have presented data on the photometric and spectroscopic variability of the moderately 
young, solar-type star k 1 Ceti. We describe our continuous photometric observations over one 
month in November 2003 using the MOST satellite and our high-resolution spectroscopic CFHT 
observations over four observing runs in 2002 and 2003. The MOST data (averaged over the 
101.4-min orbital period of the satellite) have an estimated point-to-point accuracy of ~ 0.2 milli- 
magnitude or ~ 200 parts per million, and an unprecedented temporal coverage of 30.5 days with 
a completeness of 96%. 

The MOST orbit-averaged photometry of k 1 Ceti shows two distinctly different periodicities 
- a well-defined 8.9-day variation and a smaller, less regular variation whose period appears to be 
around 9.3 days or slightly longer. The best interpretation of these variations is that two spots (or 
compact spot groups) are at different latitudes on the stellar surface, moving at angular rates which 
are different by AQ/Q ~ 4 %. This is one of the most direct evidences of surface differential rotation 
in a cool star other than the Sun and the second such after EK Dra (Messina &: Guinan 2003). All 
other measurements of SDR in other stars have been more indirect, either through season-to-season 
variations in spot period variability or through subtle broadening effects in spectral line profiles, 
so far observable only in the most rapidly rotating stars. 

The details of the smaller slower moving spot on k 1 Ceti depend on how one describes and 
removes the larger spot, 8.9 day, variability and on the treatment of the slowly varying residuals. 
The multi-parametric modeling of the spots suggests that the smaller spots is located at a higher 
latitude, which would suggest a solar-like SDR pattern. However, independent spectroscopy of the 
star, some of which was obtained only a few variation cycles before the MOST photometry, reveals 
that maxima of the Ca II K-line emission reversals are synchronized with a period of about 9.3 
days, lending support to that period identification for the second spot. While showing significant 
evolution in size, shape and/or flux contrast even during our month of observations, the 9.3-day spot 
may in fact be the more stable one and perhaps visible over decades, based on earlier observations of 
re 1 Ceti in the literature. The shorter, 8.9 day, periodicity suggests that k 1 Ceti rotates faster than 
previously supposed and thus may be even younger, based on the correlation between equatorial 
rotation period and age in solar-type stars. 

The unbinned MOST photometry, with an integration time of 40 seconds and a sampling of 
once per minute, sets limits on (1) the presence of p-mode oscillations, (2) the eruptions of flares 
over one month of monitoring of this historically active star, and (3) the occurrence of planetary 
transits in the system. The data effectively rule out p-mode oscillations larger than 25 parts per 
million (ppm) at the 3-cr limit. There were no flares with durations from 2 to 200 min whose 
white-light amplitudes could be larger than 2-3 mmag. There can be no planets with orbital 
periods shorter than 30 days larger than about 0.5 Jupiter diameter in an orbit plane which would 
lead to transits as seen from Earth. Within the value of a^y ~23 m s _1 we determine, there could 
be a close Jupiter-mass companion in an inclined orbit. 
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The information gleaned from this early set of MOST data, obtained before the satellite perfor- 
mance was optimized for normal operations on a target which was not part of the original planned 
scientific program, indicates the tremendous power of precise, nearly continuous photometry span- 
ning many weeks, which is possible from a dedicated instrument in the right orbit in space. There 
are many more applications like the ones described in this paper, and instruments like MOST can 
be an effective tool in studying the surface properties of rapidly rotating solar-type stars. 
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Fig. 1. — (upper panel) The light curve of k 1 Cet obtained from 30.5 days of continuous observations 
with MOST. The points are averages over the 101.4 min orbital period of the satellite. The 
variable background from Earthlight has been subtracted. The data are presented as relative signal 
normalized to the maximum flux observed during the run. The solid curve is the best fit of a simple 
single spot model to the dominant 8.9-day variation in the data (see Section 4.3 and Figure 4). 
(middle panel) The residual variation after subtraction of the 8.9-day model curve plotted in the 
upper panel shows clear variations whose minima are spaced by about 9.3 days (see Section 4.2). 
(bottom panel) The predicted variation of the Ca II K-line emission flux based on the spectroscopic 
observations obtained only two months before the MOST run (see Section 5.3 and Figure 9). 

Fig. 2. — (upper panel) The variable sky background measured during our observations, expressed 
as a fraction of the maximum stellar signal on JD 2,452,953. This background is primarily due 
to stray light entering the instrument focal plane from the illuminated limb of the Earth, (lower 
panels) The angular separation of the Moon from the star and the corresponding lunar phase 
during the MOST observations. There is no evident correlation with the sky background in the 
upper panel. 

Fig. 3. — The periodogram based on the original data (the points in the upper panel of Figure 1). 
The periodogram does not show any significant frequencies above one cycle per day. 

Fig. 4. — The data in the upper panel of Figure 1 are shown phased to a period of 8.90 days. The 
error bars show the ±1— <r scatter in the phase-binned data. The scatter is larger at phases when 
the second-spot variability is best noticeable, consistent with visual inspection of the light curve. 
The single spot model from the upper panel of Figure 1 is indicated by the thin solid line. 

Fig. 5. — Typical CFHT spectra of n 1 Ceti and of r Ceti expressed as normalized flux versus 
wavelength; the latter star served as the sharp line template. The broadening functions for n l Ceti 
were determined from the upper parts of the photospheric spectra after removal from the spectra 
of the lower parts close to the cores at < 1/3 of the maximum height. 

Fig. 6. — The broadening function for k 1 Ceti determined relative to r Ceti for one of the CFHT 
spectra. The broadening function gives a one-dimensional image of the rotating star, projected into 
the radial velocity space; it can be used for determination of the value of V sini. The broadening 
function is dimensionless and - for perfect match of the spectral types - its integral equals unity. 
The thin line shows the intrinsic spectral resolution of the method while the broken line gives the 
calculated rotational profile for V sini = 4.64 km s _1 and the limb darkening of u = 1.0. The spike 
in the center shows the intrinsic resolution of the technique. 

Fig. 7. — The Ca II K-line emission reversal, normalized at the ends of the spectral window in a 
single typical CFHT spectrum. 

Fig. 8. — Residuals from the Ca II K-line mean averaged normalized profile, as shown in Figure 7, 
for individual nights, with one to three observations averaged if close in time. Note that all changes 
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are confined strictly to the Ca II emission reversal. The phases of the 9.332-day periodicity are 
given in the legend; they correspond to those in Table 2 where individual observations are listed. 

Fig. 9. — {upper panel) Integrated residuals of the Ca II K-line emission expressed in A and plotted 
versus the 9.332-day rotation phase. Data points from 2002 are marked by open circles while those 
from 2003 are solid circles, {lower panel) Radial velocity deviations from the first spectrum. 

Fig. 10. — The Fourier amplitude spectrum of the unbinned MOST photometry with 40-sec integra- 
tions and one-minute sampling interval. The maximum frequency in this plot (4 mHz) corresponds 
to a period of 250 seconds; there is no noticeable change in noise level nor any signals evident 
at frequencies up to the Nyquist frequency of about 8 mHz. The amplitudes for frequencies cor- 
responding to the satellite orbital frequency and its harmonics have been removed to avoid any 
contamination by orbital modulation of the stray light background. The dashed line is the formal 
l-<7 noise level in the spectrum. 
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Table 1. k Ceti: MOST photometric data, orbital averages 



JD(hel)-2,452,948 a 


Normalized flux b 


Background flux c 


0.461914 


0.986152 


0.213726 


0.525269 


0.987092 


0.209837 


0.598022 


0.988443 


0.213698 


0.668823 


0.989125 


0.212024 


0.736816 


0.990238 


0.210084 



a The Julian Date at the effective center of a MOST orbit 
of 101 minutes. 

b The observed flux normalized to the maximum on JD 
2,452,953.0. 

c The background flux normalized to the maximum of the 
stellar signal on JD 2,452,953.0. 

Note. — The first five lines of the full table are shown here. 
The table is available only as a computer ascii file. 
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Table 2. k 1 Ceti: CFHT Radial velocities and Ca II K line residuals 



UT Date 


JD £ 




4>rot° 


Residual K c 
Flux (± 0.0017 A) 


ARVs 
km s _1 


7/28/2002 


2452484. 


.0907 


-0.110 


-0.0074 


0.0000 


7/28/2002 


2452484. 


.1062 


-0.108 


-0.0075 


0.0460 


7/30/2002 


2452486. 


.1395 


0.110 


-0.0048 


-0.0085 


8/20/2002 


2452507. 


1459 


2.360 


0.0394 


0.0104 


8/21/2002 


2452508. 


.1227 


2.465 


0.0301 


-0.0283 


8/21/2002 


2452508. 


.1371 


2.467 


0.0233 


-0.0317 


8/23/2002 


2452510. 


.1235 


2.680 


0.0085 


-0.0147 


8/23/2002 


2452510. 


.1379 


2.681 


0.0078 


-0.0256 


8/23/2002 


2452510. 


1524 


2.683 


0.0034 


-0.0084 


8/26/2002 


2452513. 


.1172 


3.000 


-0.0143 


-0.0097 


8/27/2002 


2452514. 


1330 


3.109 


-0.0222 


-0.0186 


9/09/2003 


2452892. 


.1321 


43.615 


0.0164 


0.0278 


9/09/2003 


2452892. 


.1465 


43.616 


0.0172 


0.0389 


9/12/2003 


2452895. 


.1384 


43.937 


-0.0311 


-0.0105 


9/12/2003 


2452895. 


1529 


43.939 


-0.0331 


-0.0035 



a Julian Date at mid-exposure time. Exposure times are 1200 s. 

b Rotational phases calculated using P rot = 9.332 days, T (JD he i) — 
2,452,485.117. 

Integrated residual Ca II K emission flux (in A) from a mean of all spectra 
(see Figure 8 and the text for details). 
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